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ABSTRACT 


Expert aytems are probably the rriost useful contribution of 
Artificial Intelligence. A number of expert systems are presently 
in use worldwide and many of them have come up with flying 
colors; but very little is done for developing expert systems in 
the field of numerical computations. 

Solving Partial Differential Equations (PDEs) numerically is 
indeed a complex process and needs a lot of expertise for 
choosing numerical & codiputat ional scheme, handling ei"ror control 
and stability etc. Moreover writing aird debugging special purpose 
programs in procedural languages male tlie coniplete process 
cumbersome cknd trouble shooting. 

An expert system MA, is proposed and implemented, which can solve 
quas i-1 inear , second order PDEs. MA treats the problem in a 
generic way £ind addresses the physical probleiii directly. tt is 
developed in a relations^! environment, Prolog, which does the 
problem identification, preprocessing and postprocessing work. 
Since relational languages are inefficient for performing 
procedural jobs numerical computation is done in a procedural 
module, written in C. MA directly calls the procedural module 
’intelligently*, instead of generating procedural codes. MA is 
provided with a graphics module, which gives a graphical output 
of the numerical computation. Results are stored in a database 
chain and pointers to every set of data are stored in a B+ tree 
for efficient handling of the data. MA has a knowledge 
acquisition module, which allows the user to update the knowledge 
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base 


MA provides the user a very high level of interface and 


demands minimal knowledge about the probiein in hand. An on 
help file is available for guiding users? moreover the user 
always ask why questions for clarifying his doubts. 

The implementation of the proposed system has been done 
PC/AT/XT using Turbo Prolog and Turbo C. The system operates 
MS DOS environment. 


line 

can 

on 

in 
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Chapter one ■" Itit roduct ion and Review of Object Experti'-ie 

1.1 Introduction 

With tlie increase in the computational ability of rriachii^es, the 
expectations from a user friendly scientific software is also 
increasing. The nature of the presently available software for 
computational problems has been ratlier typical - a library of 
FOF^TRAN subroutines, which follow a srictly algorthmic path and 
obviously cannot process its 'knowledge’ to solve user's problem. 

Plethora of numerical schemes along with their stability 
and error propagation properties can often confuse the user. 
Moreover partial differential equations ( PDEs ) are grouped 
into several complex h i e rar cli i cal status, each of which needs 
specific treatment to solve. A lot of oxpertase is needed to 
develop a lengthy code in procedural languages ( eg . Fo rt ran , Basic, 
Pascal etc.) for everydeiy applications in engineering desigti, 
fluid flow analysis, chemical engineering problems, heat transfer 
problems etc. Undoubtedly there is a serious need for encoding 
the PDE solving expertise. A system with a knowledge base 
containing this expertise will be able to select the most 
suitable problem solving strategy. 

A number of research groups started working with the above goal. 

In the past few years Problem Salving Environment systems for 

PDEs have been developed C1“43. They generally take the form of a 

procedural code generator along with a higher level user 
interface . 

MA is an expert system for solving PDF oriented problems. Fig.1 
depicts its overall logical structure. Present version of MA can 
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handle 'i.econd order , linear or nonlinear PDEs . Instead oF 
generating a code in procedural languages for numerical 
computation, MA 'intelligently' calls procedural number cruncfung 
modules, whenever needed. MA's conception of F^DEs and numerical 
solution methodology is realized in terms of an object oriented 
representation which is implemented in logic programming 
env i ronment . 

1.2 Brief survey of the available Problem solving environments 
for partial differential equations 

The concept of automated problem solver evolved in late 70’ s. 
Initially the research was directed towards the sole use of logic 
programming langiu^ges; but later, the inefficiency of these 
languages for procedural tasks was ejjposed. From early 80’ s 
automatic programming en v i rotmientia started coming up CI-SII. Sorue 
of the systems start with continuous description of a problem and 
others work from the PDEs which are drawn from the conservation 
laws. Finally, some systems use symbolic manipulation to derive 
all or parts of the discretization scheme while others use 
artificial intelligence techniques to classify the type of the 
probleni being solved and to fill program templates which 
implements the discretization scheme. Although a number of good 
automatic programming environments are available, the major 
drawback of most of them is the applicability for a narrow class 
of problem; this is because these systems have knowledge about 
specific kinds of problem without having the deeper 

interpretation of the knowledge. It is indeed very difficult to 
go beyond these special purpose solvers using the conventional 
techniques of software engineering. Above all, most of them 
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suffer Prom serious probienis about user interface. Coiiiplex 
fields, like numerical solution of PDEs , demand something more 
than a simple user friendly interface; user sl’iould be able to 
treat the problein and the knowledge in a symbo3ic way, which is in 
fact the paradigm of Artificial intelligence and Expert systems. 
Sytem should at first identify the problem physically and then 
only its iiiathemat i cal interpretation in terms of PDEs should be 
made. System whicli stcirt with discrete descriptions of the 
problems are likely to be more flexible but at the cost of the 

user interface, mentioned above; this also required enough user 
expertise because the user must deal with both the continuous 

problem and its discrete approximation. 


1.3 Knowledge elicitation from the world of partial 
equations 


differential 


Human mind has been searching for the truth behind the 
mystic nature from the early days of its civilization. It always 
tried to explore the nature by simulating it either theoretically 
or physically. The very birth of mathematics was endowed to the 
6'ternal search for a symbolic language, which would be able to 
describe the nature more precisely. Partial differential 
equation is a certain kind of mathematical expression, which is 
used for modeling several physical systems. 

Majority of the problems of physics and engineering fall 
naturally into one of the three physical categories. equlibrium 
problems, eigenvalue problems, and propagation problems. 
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Equilibrium problems are problems of steady state in whicl'i the 
equlibrium configuration 0 in a domain D is to be determined by 
solving tl’ie differential equation 

LC03=f — (1.3,1) 

within D, subject to certain boundary conditions 

Bir03=gi — (1.3.2) 

on the boundary of D. Very often, but not alweiys, the integration 
domain D is closed and bounded. In mathematical tenninology such 
problems are known as boundary value problems. Typical physical 
examples include steady viscous flow, steady temperature 
distributions, equilibrium stresses in elastic structures, and 
steady voltage distributions. 

Eigenvalue problems may be thought of as extensions of the 
equilibrium problems wherein critical values of certain 
parameters are to be determined in addition to the c:o r responding 
steady-state configurations. Mathematically the problem is to 
find one or more constants («), and the corresponding 
f unct i ons ( 0 ) , such that the differential equation 

LC03=o(MC03 — (1.3.3) 

is satisfied within D and the boundary conditions 

BiC0T=aEiC03 — (1.3.4) 

hold on the boundary of D. Typical physical examples include 
buckling and stability of structures, resonance in electric 
circuits and acoustics, natural frequency problems in vibrations, 
and so on. 

Are initial value problems that have an 
Propagation problems are inn,xai 
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unsteady state or transient nature. Given the initial state, the 
subsequent behavior of a system can be predicted by solving the 
differential equation 

LC0D=^f (1.3.5) 

within the domain D wiien the initial state is prescribed as 

liC03=hi (1.3.6) 

and subject to prescribed conditions 

BiC0n=gi (1.3.7) 

on the open boundaries. These problenis are also l..nown as initial 
boundary value problems. Typical physical enaniples include the 
propagation of pressure waves in a fluid, propagation of stresses 
and displaceiiients in elastic systeins, propagation of heat, and 
the development of self excited vibrations. 

Partial differential equations are classified according 

to their order, linearity, and boundary conditions. 

The order of a PDE is determined by the highest order partial 

derivative present in that equation. PDEs are categorized into 

linear, quasilinear, and nonlinear equations. Consider for 

example, the following second-order equation '■ 

a(.)dsu + SbC.)dlu_ + c(.)d£u + d(.)==0 (1.3.8) 

dy® dxdy dx^ 

If the coefficients are constants, or functions of the indepedent 
variables only C(.)5(x,y)3, then the equation is linear. If the 
coefficients are functions of the dependent variable and/or any 
of its derivatives of lower order than that of the differential 
equation C( . ) 5 Cx,y,u,du/dx,du/dy) then the equation is 
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Neumann conditions The? derivative of tlie dependent variable 
is given as a constant or as a function of the independent 
variable. For exaiciple, 

dT/d!<=“-0 at K~1 and t •;=0 

In a physical problem a neumann boundary condition can be 
specified if the gradient of the dependent variable is known, for 
e;;.ample perfect thermal insulation in case of a heat conduction 
problem may give rise to the condition dT/dx-0 at the insulated 
boundary . 

Cauchy conditions A problem which combines both Dirichlet 
and Neumann conditions is said to have Cauchy conditions. 

On the basis of their initial and boundary conditiorrs, PDEs may 
be furthar classified into initial value or boundary value 
problems. In the first case, at least one of the independent 
variables has an open region. In case of unsteady problems, time 
variable has the range 0<“t<=« and no condition is specified at 
t~« is an initial value problem. When the region is closed for 
all independent variables and conditions are specified at all 
boundaries, then the problem is of the boundary value type. 

1.3.2 Numerical solution 

The dynamics of physical systems that have more than one 
independent variable can be modelled by using PDEs. Only a few 
types of them render themselves to analytical solutions, for 
example some linear parabolic and elliptic PDEs. However, the 
majority of PDEs, especially the nonlinear ones require the 
application of numerical techniques for their solution. 
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WumeriLal solution of PDEs in itself is a wide field and a 


pletliora of schemes exists For tackling therriu No attempt would be 
made to describe all of them; only, the Finite difference scheme 
will be discussed briefly. 

1.3.3 Finite Difference Approximations 

The calculus of fiiiite differences may be characterised 
as a "two-way street" which enables the user to take a 
differential equation and integrate it numerically by calculating 
the values of the function at a di s c r ete ( f in i tc > number of points. 
Iti differential calculus, the definition of the derivative is 
given as 

df(i<)/dK 1 j.;o =f ' (!;o)"lt f (!C)-f (xo) -"( 1.3.10) 

X~>XO X “ xo 

In the calculas of finite di f f e r en ces , the value of (x-xo) does not 
approcich nero but remains a finite quantity. If we represent this 
quantity by h, 

h = X -■ xo 

then the derivative may be approximated by 

f • ( xo >5 f (x)-f (xo) (1.3.11) 

h 

This concept of replacing the continuous operators by discrete 
ones can be extended for backward, forward and central difference 
approxiniat ions by using Taylor series as shown below: 

Development of Taylor series for u(x+5x,y) about (x,y) gives 

u(x+<Sx,y) = u(x,y) + Si< ^(v y) + (<5x^/2! ) dgu f y) 

dx " dx2 

+ ( (5x)^/3! )d~U(;<^yj + 0<<5x'^) 
dx^ 
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which upon division by fix, results in the relation 

iLii(x,y) " ^.u(x-^(?x,y) ■■ u ( x , y ) IJ/dJx + 0(<5x) - - f 1 . 3 . 12 ) 
dx 

This is known as the forward difference appro irnat i on of the first 
order differential operator. Similarly using the Taylor series 
for u<x--<Sx,y) about <x,y) backward and central difference schemes 
can bo obtained. 

Table 1 shows a listing of the Fijute di ff erence approximations 
of partial derivatives using central di Ff er enoes . 

1.3.4 Stability and convergence 

The presence of round-off error or any othc'r 
computational error may lead to numerical instability. 
Consider the simple ordinary differential equation 

dy/dx - y - X — (1,3.13) 
whose general solution is 

y =- Ae'^ + (x+1) —(1.3.14) 

where A is the integration constant. In an approimate solution 
the eponential term is likely to be introduced as a result of 
round-off errors. Any numerical scheme which allows tlie growth of 
error, eventually 'swamping ' the true solution, is unstable. 
These numerical phenomena must be avoided by restrictive action, 
such as limiting the interval size or adopting an alternative 
method , 

Let us talk about stability in a more mathematical way. Let 
U(x,t) be the solution of a given difference appro imat ion , 
solvable step by step in the t direction. The effect of a 
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inistake or, more likely, round off error in machine computation 
may replace UCiCort^) by UCKort^) + e at the grid point (KQjt^). If 
the solution procedure is continued witii the value U(;co,tQ> + e 
without new errors being introduced, and if at subsequent points 
the value U (;;,t)is obtained, then we denote by U ( ;< , t ) -U ( ;< , t ) 
the departure of the solution resulting from the error e at 
(Ko,tol- When errors are introduced at more than one point, 
cumulative departures result which are not additive ecept in 
linear problems. If 6 is the maximum absolute error — and h 
thie interval size, tlien a procedure is said to be pointwise 
stable if the cumulative departure tends to zero as S — .-=0 and 
does not increase faster than some power of h ^ as h — >0. 
when the corresponding continuous solution remains bounded a 
finite difference process within the semi"-inf in ite strip 0<x<1, 
t>0 is called stepwise unstable if for a fied network and fied 
homogeneous boundary conditions there exists initial disturbvinces 
for which the finite difference solutions become unbounded as 
j — >« . 

A second fundamental concept, that of convergence, is often 
related to stability. To introduce this idea we utilie the PDE 
L<u)=0 in a region D, u=g on T 

where T is the boundary of D- Let us suppose that there is 
only once, say h, and write the finite difference problem as 
Lf-, (U)=0 in D, U=gi-, on T. 

We say the finite difference scheme converges if U(P) 
converges to the solution u(P), with the same boundary values, as 
h — >0. 
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1-3.5 Integration in time domain 

Integration of the dependent variable can be approximated by 
a weighted bacl::ward difference along the time axis. Depending 
upon the weighting factor several schnies are proposed. 

For example, a second-order j.>artial derivative can be expressed 
as a weighted average of tfie centraxl differences at points 
and (i,j), where a denotes the spatial directiotj and j 
denotes tiriie di recti oaa s 

d®u/dx® — + ( 1 ""G) !?^u ^ j 

= ec ( 144..1 ^ - 2ui^jH..| -r Ui- -I ^ j 4 .^ ) /h2;i + 

+ j - Sui^j + ^ j )/h2D 

- - (1.3.15) 

where 0 is in the range 0<,9:£1. Sindlarly other differential 
operators can also be integrated with a time weighting factor. 

Explicit Scheme s 

When ©=0, equation 1,1 reduces to its second part, 
which involves property values at the previous time step. This is 
known as explicit scheme. The stability of this scheme is 
conditional . 

Crank-Nicolson Scheme • 

When 6=1/S, equation 1.1 yields the Crank-Nicolson 
formula for second order differential operator. This is also 
conditionally stable. 

Implicit Scheme • 

Finally, when ©=1, equation 1.1 becomes the implicit 
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f orniulat ion . Implicit formulas are uncoiidi t ional ly stcible. 

It can La generalized tliat most explicit finite difference 
aiJproKimations cxre conditionally stable whereais most implicit 
apprvjK imat i ons arc* unconditionally stable. The explicit methods, 
however, are computationally easier to solve than the implicit 
techniques . 


1.4 Present work 

Several limitations of the present problem solving environment 
(PSE) for partial differential equations, as mentioned in section 
1.2 iiave motivated many researchers C1-6S to work for systems 
witli greater flexibility for handling equations along with a 
relational handling of the vast knowledge base. 

We propose a system, liA which has several unique features for 
solving PDEs and equipped with some powerful AI techniques. MA 
eploits the properties of both relational and procedural 
languages. It has a robust PDE solver along with a higher level 
user interface. F i g . 1 . oJofH shows a schematic diagram of the 
system. The main goal of the expert systems is to treat the 
problem in a natural symbolic way. MA aspires this by handling 
the PDE solving problem as a physical problem and is not limited 
to the computational procedures. 

MA consists of three main modules t 

1) Preprocessor 

2) Solver 

3) Postprocessor 

The logical organization of MA is shown in Fig.1. 
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Preprocessor modulp identifies the problem with the help of a 
physical problem knowledge base and receives required parameters 
for defjning it uniquely. It also points out errors in typjng. 
User defijied paraineters are checked for stability criterias, 
specified in the rule base. 

Solver is a procedural module which does the numerical 
coinputat i on for solving tl’ie set of algebraic equations, generated 
from the discreti::ed PDE, which is written in C. This module 
directly does the computation work by calling suitable number 
crunthiJig procedures, instead of generating a code in procedural 
language and then mal.:ing it an executable file, which is a 
conventional feature for most of the automatic PDE solver. This 
unique feature makes MA more competent for microcomputer usage. 
Postprocessor module handles the solution kept in database chain. 
In case of time dependent problems times are kept in a B+ tree 
chain along with a pointer to the corresponding solution values. 
This makes the search for the solution at a specific time from a 
large set of solutions easier and fast. Adequate graphics 
facilities are available in the Postprocessor for plotting 
different kinds of graphs - 

MA treats the problem in a generic way and is not 
restricted by any ad hoc implementation strategy. This will be 
described in chapter 3 in details. 


20 



Cl)aptEr ClMO 


Logic Programming 

Prolog & 
Expert Systems 





Chapter Two - Logic Pr og raiiuning Prolog & Eiipert Systems 

2.1 Introduction 

Logit pr ograriiming is an approach to computer science in 
wh;( ch the horti clause form of the first order predicate logic is 
used as a high level prograrnniing language. The study of symbolic 
logic goes back to the work of Aristotle in the fourth century 
B.C. First order predicate logic is a branch of symbolic logic 
that has largely evolved in ttie twentjetli century. Logic 
programming is based on a subset of first order predicate logic, 
but IS equally broad in scope. Logic programming allows a 
programmer to describe a situation with formulas of predicate 
logic and then to use a mechanical problem solver to make 
inferences from the formulas. 

2.E History of symbolic logic 

What is now known as "traditional logic" began at the time 
of Aristotle over 22 centuries ago. Aristotle attempted to codiPy 
into a scientific system the way that knowledge can be most 
effectively pursued through rational debate. 

Aristotle analysed the form of a statement into the following 
e 1 ements s 

Quantifier Subject Copula Predicate 

so that, for instance, in the following statement! 

Some IITK students are intelligent 

"IITK students" is the subject, "Some" quantifies the subject, 
"are"is the copula and "intelligent" is the predicate. 



Traditiona] logic recognises -Tour possible relationship between 
two classes. They are: 

All A is B. 

Nu A is B. 

Some A is B. 

Borne A is not B« 

where A stands Tor a subject class of items and B staTids for a 
predicate class oP items, llel 1 en ist i c Greece also talked about 
syllogisiri. Tite syllogiHin is a set of rules which govern the 
(.onclusion from a set of preiruses. But Aristotelian logic was 
incomplete and inadequate to describe real world situatiojis. It 
was the English mathematician Deliorgan and Boole who criticised 
Aristotelian logic first in the middle of the iiineteentli century. 
Although traditional logic is a logic of classes, it contains no 
notion of complement. Moreover its complete dependence on natural 
language has led to many absurd difficulties. 

Boole and DeMorgan associated logic with mathematics rather than 
pliilosophy. They developed the concept of set theoretic operators 
for logic. That was the birth of modern symbolic logic which led 
to Propositional and Predicate logic. 

E.3 Propositional and Predicate logic 

2-3.1 Atomic 

Formulas 

A proposition is any statement that 
can be assigned a truth value (either true or false). In the 
conventional lexicon of propositional logic, atomic propositions 
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are defined 


propoHitians that cannot bo broher, down into 
components. In F’ropos i t ional logic, tl.e basic truth-valued object 
ill an atomic proposition that cai\not bo subdivided into 
coinponenls in any way. Tin; basic truth valued object in predicate 
logic is an cxlomic formula. An atoriiic formula is composed of a 
predicate symbol, and termis that act as arguments to the 
predicate symbol. In general, the predicate symbol Is the liame of 
a relationship that holds between the arguments. As an atomic 
formula is itsulf constructed of otl'ier objects, it posseses mucii 
greater expressive power than an atomic proposition in 
propositional logic. 

An atomic formula is written v-^s a. predicate symbol followed fay 
some number of arguinents inside parentheses. Each argument is a 
term. The general form of an atomic forniula is s 
P(t1,t2 tn) 

where P is a predicate symbol and t 1 , tS tn are terms. 

Term! 

A term can be one of the three things : 

1) a constant? 

2) a variable? or 

3) the application of a function. 

The application of a function is written as a function symbol 
followed by at list of arguments inside parentheses .Each argument 
is itself a term. The general form of the application of a 
function is s 

f (t1,t2, tn) 

where f is a function symbol and t1,t2....tn are terms. 


23 



2.3.2 Well formed 

formulas 

The result of conibiniug atomic formula together with 
logical connectives is a well-formed formula (wff). A well formed 
formula can be defined as follows s 

1) an atomic formula is a wff? 

2) if A and B are wffs then so are 

Wff Read as 

"^A not A 


A & B 


A and B 


A V B 


A 0 r B 


A — B 


A implies B 


A • ; B 


A if and only if B 


For some 


Vx A 


for all 


Quant i f i e rs 


and V are known as existential quantifier. 


As in propositional logic, the truth value of a wff 
depends on the interpretation .An interpretation of a wff can only 
be made with respect to particular domain of interpretation, 
which is the set of all possible values of terms that occur in 
the wff. 

2-3.3 Inferences in 
predicate calculus 

Inferences in the propositional calculus result from 
the internal structure of sentences that are comprised of 
component sentences. However, the inference in predicate calculus 
needs an analysis involving a finer breakdown of structure than 
one stopping at atomic sentences. Let us consider an example for 
a clear understanding of the subject-predicate relationship. 
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varibl^ instantiation and existential quantifier. 

All kids are nice. 

A-1 IS a kid. 

Hence, A‘1 is nice. 

Let kids(y) indicate the fact that the object y Is a kid, and let 
niceCx) indicate the property that x is nice. Both kids< ) and 
nic.e( ) are called predicate of arity 1, since each involves one 
c^r gudient . 

using the notions thus far introduced, the inference (1) tan be 
byiTibo 1 i zed by 

V K C k i ds ( X ) — > n i c o ( x > n 
kids ( A‘1 ) 

Hence, nicr»(A1> 

Grammatically, the instantiated predicate nice(A1) represents the 
declarative statement "A1 is a kid." and is either true or false, 
verity is based on the distribution of facts. But what about the 
expression kids(x)? Factually, the expression is empty and is 
not a proposition. The set of objects that the variable in an 
open sentence might represent is called universe of discourse. 

E.4 Prolog - PROgramming in LOGic 

After discovering the computability of logical 
consequence in an abstract sense, it is a natural step to want to 
mechanize the proof process on a computer. Although full 
predicate logic is a very expressive language, implementation of 
proof procedures of predicate calculus in a general way may cause 
combinato r ial explosion. In 1965 Robinson came up with 
resolution, an inference rule appropriate for machine inference. 
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Loveland, Kowalslci, and Kuehner -further refined this wath the 


techniques of model elimination and tl'ie seiection function. In 
early -^-eventies. Col rnerciuer and Roussel brought this work 
together into a language that they called PFJOLOG (for RROg ranicning 
in LGGi r ) . 

Following parts of this section describes resolution 
technique, which is the logical basis of the language, and 
severcil other features. 

S.4.2 Resolution 

Resolution worts as follows! 

two clauses can be resolved with one another if one of them 
contains a positive literal and the other contains a 
cor responding negative literal with the saxme predicate synibol and 
the same number of arguments, and if the arguments of both 
literal can be unified (matched) with one another. The clausal 
forni of predicate logic is a way of writing formulas that uses 
only the connectives &, v, and A literal is either a positive 
or negative fornriula. Consider a theory composed of the following 
two clauses! 

P(a) V '^Q(b, c) (1) 

Q(b,c) V '''R(b,c) <E) 

Since clause (1) contains the negative literal '"Q(b, c), and 

clause (2) contains a corresponding positive literal Q{b,c), and 
since the arguments of the two literals can be unified, then 
clause (1) can be resolved with clause (2). The resolvent is 
shown below ! 

P(a) V '^R(b, c) (3) 
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Now, aJ.1 the three clauses can be used in future resolution. 

T!ie basic problem is to prove whether a clause is or is not 
a consequence of a theory. To be able to automate the discovery 
process, we would lilce to find an efficient rule tliat can detect 
the inconsist€;ncy of a set of clauses. 

Tliere are more than one problem solving strategy that can be 
pursued using the resolution rule. Prolog incorporates a top-down 
(or backward ) strategy. Tliis strategy aims to detect whether a 
single clause, C, is a consequence of an existing set of clauses, 
T.The set of clauses, T is assumed to be consistent. The 
algorithm works as followss 

To begin with the negation of the clause is being tested, ~C , and 
is added to the existing set of clauses, to form a new set of 
clauses, T’.If the algorithm can derive the empty clause from T' 
it is inconsistent because of the presence of '^C , and C must 
therfore be a consequence of T. 

2.4.3 Prolog syntax 

The Prolog language is a combination of powerful ideas, 
including '• 

1) the use of Horn clauses to represent knowledge? 

2) descriptive style of programming? 

3) both declarative and procedural semantics. 



Ill fA Horn clause, one conclusion is followed by zero or more 

condi lions, wliich is wrilten in Prolog as Poliows = 

conclusion s ■" 

condition 1 , 

conditions, . 

conditionN. 

is read as "if", and , is read as "and"? so the whole clause 

can be read s 

The conclusion is true j f 

conditioni and conditioiiS and. . » cond i t i onM 
are all true. 


2.4.4 Semantics of Prolog 

The semantics (meaning) of a formula of a symbolic 
logic. refers to its truth value. The semantics of a constant 
symbol inside of a formula refers to its value with respect to a 
domaiin of interpretation. In computing, on the other hand, the 
semantics of a programming language construct usually refers to 
the behavior of the computer when the construct is evaluated. 
Because Prolog is both logic and a pr ogramming language, both 
notions of semantics are applicable. 

There are three seniantics niodels to explain the meaning 
of a Prolog program, which were first explained by Kowalski. They 
are ! the declarative model, the procedural model, and the 
abstract machine model. 

The declarative semantic model of Prolog specifies the truth 
value of the relations. The word declarative is used because a 
Prolog clause declares that a relation holds between its 
arguments if all of the conditions of the clause are met. 
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Fof instance, the foilowiiic) clause s 


s:<et.utiveCNMe, Salary) !~ 

employee (Name, Salary ) , Salary > 3000. 

can bo read as i 

Atiyone (Name ) is an exetutive if 
I'lO or she is an employee 
with a salary gre^itor than Rs 3000. 

Read according to tlie declarative model, Prolog clause are 

formulas of first order predicate logac written iri Horn clause 

form. 

According to the procedural model, tiie conditions of a clause 
specify a process to establish a truth value of the conclusion of 
the clause. A set of clauses wjth the sairne predicate nanie and the 
same number of arguments as a proct'dure is understood to be a 
call to that procedure. Recid procedural ly, the meaning of the 
above clause is : 

One way to find an executive is : 
first, find an employee, 
then second, verify that the salary 
of the employee is greater than Rs.3000. 

A Prolog query has a computational meaning in the sense that they 
trigger a certain behavior on the part of the Prolog interpreter. 
The interpreter applies a problem solving strategy to evaluate a 
query against a set of clauses? its problem solving strategy can 
be characterized computationally by an abstract machine. The 
abstract machine model specifies the meaning of a query and set 
of clauses in terms of abstract machine actions. 
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2.5 Knowledge representation 

One of the main goals of the AI rei>eetrch coinniun i ty to js 
develop system which can simulate the real world knowledge in a 
machine. The art of encoding ’knowledge’ is technically !:nown as 
l:;nowledge representation. Wo have already discussed the 
representational formalism of first order logic. Now we shall 
discuss several means of representing knowledge briefly. 

2.5.1 Facts and Rules 

As we mentioned earlier, the a Prolog program consists of a set 
of clauses, which is either a fact or a rule- A fact is an 
assertion that a particular relation holds; it is written as a 
name followed by the arguments inside a piar entfies i s . The 

following fact expresses the idea that "ma solves pde" s 

sol ves (nia, pde ) . 

A rule is a fact whose truth value depends on the truth values of 
other facts. A rule can be represented in a Horn clause with each 
clause using assertions as well as the rules in the knowledge 
base. An example in Prolog is as follows s 

X Y , Z. If Y and Z , then X. 

X Y ; Z. If Y or Z , then X. 

2.5.2 Semantic network 

In a semantic network, entities and classes of entities are 
identified with nodes, and relations between entities are 
identified with arcs joining nodes. An arc connected to a single 
node establishes a property of that node. Several types of 
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attributes are often employed, for instance the clcxss relation 
:is-v^ and the possesation relation has-a. In 0-A--V 
< obj e c t-"at t r 1 hut e-"val ue ) notation examples ofthis relations may 
loo k 1 ike s 


i s__a C pde , equat i on J . 

is_a(ode, equation) . 

has__a( pde , two_inde pendent_var iables ) . 

An essential concept of the semantic network formalism is that of 
hierarchy, which makes it particularly adept at representing 
taxonomies of knowledge. Each level of the taxonomy is 
represented by a node connected by is_a arcs to tiigher and lower 
levels . 

2.S-3 Frames 

Minsky L;133 originated the concept of a frame as a way to represent 
knowledge about situations. Each frame contains slots thoit 
identify the type of situation or specify the parameters of a 
particular situation. Moreover, other forms of knowledge, such as 
procedural information, can be used in slots. Inheritance is one 
of the important Features of frame formalism. It is possible to 
specify that if a slot in one frame is not specifically filled, 
the frame will inherit a default value of that slot from a 
superior frame. Inheritance of slot values is enabled between two 
frames by the presence of an a_.kind_of slot in one frame that is 
filled by the name of the other frame. The a_kind__of slot is 
similar to the is a arc in a semantic network. 
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E.6 Expert systems 


Expert systems are a new class o-P application programs, inteiidod 
to ifial-'e ti-ie knowledge op an expert in some special field readily 
available to ''lay'' users via a computer driven, interactive, 
dialog oriented systeiri. The fundamental differevice with respect 
to traditioTial pr og ramming is that the processing rules are 
stored in a knowledge base along with the data itself. Tho cot'e 
of £in expert systeni is an inference engine. It applies a 
particular strategy to draw conclusion from the knowledge stored, 
theroby producing new knowledge. Suitable heuristics must provide 
for targeted knowledge processing, so that a meaningful response 
to cl user query can be produced within a reasonable amount of 
time. 

Fig. 2 shows a schematic diagram of a general expert system. 
Modern expert systems comprises two other important features! 

1. knowledge acquisition module 

2. an explanatory interface 

The knowledge acquisition module can be seen as a kind of back- 
end to the knowledge base; and the explanatory interface as a 
front-end to the inference engine. But perhaps the most fruitful 
way of looking at this diagram is as a pointer to the four 
central unsolved problems in expert systems methodology. 

These are! 

1. the problem of knowledge representation 

2. the problem of approximate reasoning 

3. the problem of knowledge acquisition 

4. the problem of human-machine interfacing. 
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It may «ieHni parado:; 3 c.al to iay thait tlie "itatw of art in oKpsrt 
systems is a collection of unsolved probieiris, when hundreds of 
firactical expert, systems are in worldwide, but it is true. 
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Cl'i-Hpter Three - Gy stein Analysis 


3.1 Introduction 

For tiie numerical solution of PDFs there is the never 
endin '3 discussion of special purpose program versus general 
"blacl- boK" solver. An ideal Dlauk bo:i solver should have the 
fallowing properties ! Full flexibility of 

1) PDF. operator < all types of linear nonlinear systems of 
PDFs) , 

P.) boundary condition operator ( all types of nonlinear DCs), 

3) solution operator < variable quality, error estimate, 
se 1 f adaptat ion ) , 

4) geometry ( arbitrcxry complicated 2~D or 3-D domain ). 

As we discussed in cliapter 1.2 all available solvers are 
compromises. Our black box solver MA has the following properties 

1) solves quasilinear systems of PDFs, 

2) admits arbitrary nonlinear DCs, 

3) restricted to a rectangular doinain <or a domain that can 
be transformed analytically to a rectangular domain ). 

Full geometric flexibility can only be obtained by FPIM on 
unstructured grids, but then we have a data structure that 
ne ccess i tates much indirect addressing and thus reduces the 
efficiency of vectorization. Above all we have for the FEM, 
problems with the requirements 1 to 3 mentioned above. The FEM 
thus needs always an individual adaptation to a distinct problem. 
This was the main reason for choosing FDM. 
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3.2 DOMAIN OF EXPERTISE 


Expert system MA was initially designed to serve as the 
computat i onal front end of a knowledge based system for heat 
exchanger design, Even then, MA is flexible enough to treat its 
problem domain oF solving PDEs numerically in a generic way to 
other applications. 

At present it can handle second order PDEs. This covers a 

wide range of CMigineering and scientific appl i cat i ons . These 

equations can i^e ex. pressed in the fonn - 

A*d?u/dx^ + B«dSu/dxdy i- C^id^u/dyS + D*du/dx + E*du/dy 
i- F =••• 0 “(3.2,1) 

where A,B,C,D,E,F and const are either constants or functions of 
any combination of x,y,u; u is the dependent variable and x,y are 
independent variables. Non-linear terms are linearized by 
Jacobi's method using Taylor series. 

MA must be supplied with the required boundary conditions and 

possibly initial conditions in order to obtain a unique solution. 

It can hc^ndle boundary conditions in the form of Dirichlet , 

Neumann , Periodic or mixed conditions. After recogi\izing the 

problem with the help of its knowledge base it chooses an 

appropriate numericai algorithm and then goes for the solution. 

Knowledge base of MA contains facts and rules which it uses to 

achieve goals iir the process of solving PDEs. It can always be 

updated and an experienced user can control the algorithmic path. 

Knowledge base rules contain the fallowing : 

13 Classification of PDEs 
23 Convergence properties 

33 Applicability of discretization schemes 
43 Choice of computing scheme 
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SuiTiG of the fact:; are stored in frames . Some of them are s 

13 PDEs and their corresponding problem domains 
2.2 Di s c r et i ;.':at i on scheme:^ 

3TJ Grid information 

43 Inf ormatiovi about coniputat a onal s<-hemes 

3.3 KNOWLEDGE REPRESENTATION 

In order to reason about the algoritl'im chosen to solve the PDE , 
MA irtus t have a knowledge base cotitaining the expertise on problem 
domiUJi . MA has an object oriented knowledge base about the world 
of PDFs , stored in rules and frames . We shall discuss the 
knowledge representation in detail here- 

Prolog implements a subset of first-order logic , the Morn 
clauses . A Prolog clause has both procedural and declarative 
interpretation . Constants , variables and structures can be used 
as argurneiits to predicates . More than one Horn clause may be 
needed to define a predicate and a cor respoirding AND/OR graph is 
thus created for each predicate. 

3-3.1 Rules 

The main clause of MA which controls the process of solving PDEs 
right from recognizing it can be stated as s 

s tar t._pro cess if 

get_problem(Equat ion , Type), 

get_gr id__ parameters CTy pe , Param_l i st ) , 

discretizeC Equation ,Param_l i st , Dis c r ete_l ist ) , 

SQlve(Discrete_list) . 

start__proc8ss is a dummy predicate which fires a chain of actions 
, finally solving the PDE ; get_profalem receives the problem from 
the user and then instantiates variables - Equation ( a list of 
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'strings ) and Type (a structure containing the identifiets oT the 
eguation). These inherited terms are passed in 
g f* t^g r i d„pa r amt.; t e r s , d i s t.. i- e t i e and solve. 

3.3.2 Frames 

Frames are used in MA for storing st raigi'itf orward informations. 
A frame is composed of slots for variables and their values. Here 
IS a frame for storing discretization schemes - 


D i s c !" e t i z a tion scl'ieme 


nch€?nie1 I 

A 

t 

1 

schemeH I 

B 

t 

‘ . « * 

1 " 1 

* • 

« t 


The corresponding representevtion in Prolog in given below t 

frame ( di s c rot i zat i Qn_s theme , scheme 1 ,a) . 
f re^iiie < di s c r et i zat i on_s c heme , s chemeE, b) . 

Here s cheme 1 , s chemeE . . represent the names for Uie slots such as 
implicit , f ul 1 y_e:<pl i cit etc. and a,b.. are structures containing 
cor res ponding appro:<imat6d(discrete) terms. 

3.3.3 Object oriented representation 

Object oriented hierarchial represention has many remarkable 
advantages in contrast to the single stage classifier systems. MA 
classifies tl'ie problem from three aspects -- mathematical , 
discretization and computational . 

Fig, 3 depicts the hierarchical relationship among the objects 
based on mathematical taxonomy. A sequential search along the 
three main branches ~ equation , spatial domain and time domain 
finally identifies the object . Similar hierarchical relationship 
exists for the discretization and computational schemes. 
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3.4 USER INTERFACE 


Usor interface is probably the most important factor for the 
survival of a ■..>oftwai'e . MA has an interface designed to interact 
at a higher level tlirough pull-downs , pop-up-metjues , editors 
and at the same time to be transparent enougli to have a ’feeling* 
about the ongoing internal process . 

Fig. 4 depicts the main menu of MA. Two options are available 
for etitering the problem — 

Feed equation s ■■ User can directly feed the PDE to be solved 
through an editor. 

State problem Instead of explicitly mentioniing the equation to 

be solved the continuous physics problem can also be stated , 
which in turn gathers necessary informations from the knowledge 
base . 

Default denotations of the dependent and independent variables 
can be changed by using option menu . 

Any error in typing the equation may cause unwanted errors in the 
computation. In order to prevent this , MA is provided with a 
smart parser which can point out the mistake. At any stage 
user can interrupt the system and enquire about the process. 

Plethora of discretization and computational schemes often 
confuse an .inexperienced user. MA always offer the best scheme to 
solve the PDE , but the user can also get rid of this. 

Graphics is an essential feature of MA. Its graphics modules has 
several features. 

Graphical output s 

Graphical representation of the result is often more 



c oniprohenn i VR than a simple set of numbers. MA is provided with a 
yraph plotting iiiodule which can also be called sc^paratel y . Fig. 5 
shows ooiue graphical output of the system. Usual features for 
providing legends and labels are available. At present the 
coniplete graphics module is written in Prolog. Handling and 
processing of solutions are unique and designed carefully for 
accessing tlierii efficiently. All the solutions are kept in tlie 
foriri of ©eternal database, alongwith an index kept sorted in B+ 
trees with a pointer to the cor respondiJig data. 

3.5 Controlling inference s- 

Bacl:.wi»rd tliaining is one among the primary means of performing 
inference. In this scheme verity of a syllogism is checked by 
making certain that in each of the instances in which the 
conclusion is false at least one of the premises is also false. 
This is par t i culsir ly suitable for our purpose , ais the number of 
subgoals at every level of the hierarchical tree is limited. This 
process continues until either the goal is proven or failure 
occurs , in which case system backtracks to find alternate 
solutions. Object oriented hierarchy mainly controls the search 
path. All the properties of the superclasses are inherited. Aiter 
defining the controlling path fay traversing along the 
hierarchical tree, delicate control is performed by rnetalogical 
control. A sophisticated scheme for controlling the flow path at 
run time is implemented , which makes use of dynamic cut , 
assertion and retraction of informations about the problem such 
as type of equation , grid parameters etc. 
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3.6 Equation approximation and computation 


MA dete rrnine'j a] 1 the chat' acter isti cs of an equation shown in 
Fi'j.'l by iiieanu of passing the equation recursively down through 
tlie three main branches of tl'ie hierarchy. Fig. 6 shows -soitie of tiie 
f ei'd-backn to the user and the dis(:retir:ation menu. liA solves 
PDFs by finite difference methods. For transforming the 
continuous domain PDF into a set of algebraic equations in a 
discrete domain , differential operators ritust he approximated 
MA of Per "i several schemes for discretizing the equation , eg. 
CrancF-Ni Colson , F'ully explicit etc. for time integration and 
use central diff rente approx iiuat i on for the spatial derivative 
terms. MA always checks for the convergency criteria of the 
chosen scheiiie against the grid parameters. It warns the user 
about the possible trouble and offers a chance to modify the 
parameters. If the user defines his own scheme by updating the 
knowledge base, he can always define the stc^bility criteria for 
choosing appropriate grid parameters. 

Logic: programming environments are not that efficient to serve 

procedural jobs. MA incorporates a set of procedural functions 
written in C. Linear equations can be handled in a straight 
forward way; but for non-linear terms C-modules in turn call a 
symbolic coniputat ional module to simplify the terms. 
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Clmi.'it'i' Four Imp I enieulat L on 

4.1 Int roduct ion 

So far we have discusned the baeit. knowledge of nutiierical 
aoJutioji of PDEs, concepts of Logic programming and a conceptual 
analysis of tlie proposed system, MA. In this section we shall 
talk about the practical implementation tecl'miques and some 
spi‘C3<il features of MA. 

4.P. Why Prolog ? 

We have alre<'Hiy seen some of the powerful features of 
Prolog J.n chapter II. Although Prolog is certainly not a very 
good language for procedural tasks such as solving simultaneous 
algebraic. equations, matrix operations etc., it can be 
effectively used to develop a higher level interface for a 'lay* 
user. In fact it is the preprocessing part which needs an 
relational environment for identification and initial processing 
of the problein. Prolog is one among the best suited languages for 
this kind of job. We exploited the inherent backward chaining 
feature, backtracking and pattern tnatching to develop MA. 

Major problem with procedural languages is that the mecxning of a 
construct is confined to itself and it is very difficult to 
develop a relational formalism among them, which is a inust for 
identifying the problem in a hierarchical search space. 

4.3 Why Turbo Prolog ? 

Turbo Prolog is a comperatively younger version of 
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P 1 O j. O 9 . 


l-.VQll 


irujii , 


cl nuifiber of useful •Tfcitur e s h.-ive iii-^ide it 


attract 1 ve fur developing expert systems- Followings are some of 
the importcint character isti c s of Turbo F’rolog » 

1. Turbo Prolog is descriptive 

Insttsad of a series of steps specifying how the 
computer inuot work to solve a problem, a Turbo Prolog 
programni consists of a description of the problem 
data. 

/i. Foreign language calls 

Every loinguage has its own limitations. Turbo Prolog 
is also not an exception. FORTRAN is good for 
numerical computation, C and Assembly are handy for 
for lower level access of the machine. All this 
resources can be exploited by using Turbo Prolog's 
interface to C, Assembly, FORTRAN and PASCAL. 
Similarly Turbo Prolog routines can also be called 
from foreign languages. 

3. Graphics 

Turbo Prolog supports the full featured Borland 
Graphics, which has a strengtti of 60 different 
graphics predicates. 

4. Internal and external database 

Turbo Prolog's internal database can be used for 
storing informations in RAM at run time, which is 
limited by the hardware, usually upto 6'^IOK bytes. 
Externa] database is another feature which makes 
Turbo Prolog efficient for handling large amount of 
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data. Go I'lh i fit i «. ated acl'ieirieii are available -for storing 
d a t a i r i a b i nary t r e e . 


4.4 Parsing 

F^arsing in an important part of MA. Since user is 
allowed to present the probleni syiribo J i cal ly a robust parser is 
needed to identify t!ie problem. Once the string input is received 
a binary tree representing the equation is generated. A preorder 
search acomplishes the equation identification job by pattern 
matching witl'i pre-dofined templates. Terms are simplified by 
calling a predicate simp of arity 2. Tlie simp clause is shown 
!■) ere : 

simp( lnitit\l_l ist , Final„list) !- 

1 ist^exp( Initial_list, Init ial__ei<p) , 
reduce ( In it ial_ex p , Fi nal_e:; p ) , 
e ;< p_l i s t ( F i na l_e ;< p , F i na 1_1 i s t ) . 

Here predicate list_eKp converts the string list to a recusive 
bitsary tree Initial_e:<p; reduce simplifies the binary tree 
Jnitial_eKp to Final_e::p and eKp_list converts a binary tree to a 
string list. 

Let us cosider an example! 

Expression to be parsed s (2»uim1-4-»ui+2*uip1 > /Dx 

At first this string is converted into a list of tokens. List__e;<p 
creates the binary parse tree of the list, which looks like. 

di V ( plus (minus (mult( real (2> ,var{"uim1'’) ) , 

mult C real (4> , war ( "ui") ) ,mult ( real (2> , war ( "uipl") 1 , war ( '‘0x**> > 
Predicate reduce, which does the symbolic computation, will be 
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di 5 cu<b&ed in the following section. 


4.5 Symbolic Computation 

Byml;olic: computation is necessaiy as the form of the 

given equation is unknown and the functions are of unimown 
c oriib i licit i oil uf the primitives. Symbolic coniputation is done in 
a recursive pr o cedure , which handle strings in tlie form of binary 
tree. All thn mathematical functions sin, cos, tan, cosec, sec, 
cot. In, exp are available for use. MA has two symbolic 
coriiiHitat i on modules one is written in Prolog and the other is in 
C. Ihe main problem of calling procedures written in Prolog froni 
C is that ttie heap tiiemory used by these procedures is not 
released even after they succeed or fail. This is because, the 
C iriodule controls the memory management while calling Prolog 
proteduies. As a result the system soon runs out of memory wlii ch 
puts a severe limitation on rigorous number crunching. This 
problem is solved by developing a symbolic computation module in 
C, which liandles a similar data structure. In case of nonlinear 
coefficients this module comes in much of use. 

We shall explore the reduce predicate, mentioned in the above 
section, now; reduce simlifies the expression by a preorder 
traversal of the parse tree. Look at the following primitive 
symbolic computation module, written in Prolog, which allows 
operators ", "/" only. 

/# Addition */ reduce ( plus ( X , Y) ,R > !- 

reduce<X,R1 ) , 
r6duce<Y,R2) , 
plur(R1,RS,R) . 
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/■!! subs tract ion ■»/ 


/ Si nm ] t i p I i cat i on */ 


/ -Si- d i V i i> .1 0 n 5t' / 


/s! default •!!/ 


reduce (minus < X , Y) ,R ) s™ 
reduce ( X,R1 ) , 
r educe <Y,R2) , 
rninr(R1,R2,R) . 


reduce<iiiult(X,Y) ,R) s- 
reduce C X , R 1 ) , 
reduc;e(Y,R2) , 
mulr(R1,R2,R) . 

r educe ( div ( X , Y) , R ) s- 
reduce ( X,R1 ) , 
i'educe(Y,R2) , 
divr(R1 ,R2,R) . 

reduce (R ,R ) « 


plur C real < X ), real (Y) ,R ) 

R•~X^-Y. 

rrdnr ( real ( X ), real (Y ), R ) 

R=X--Y. 

mul r { real ( X ), real (Y) ,R ) s- 
R-=X*Y. 

di vr ( real < X ), real (Y) ,R ) 

R=X/Y. 

reduce calls itself recursively till its arguments become the 
pr iiiii 1 1 ve ,» var < . ) or real<.). Then plur, minr , mulr, dive do the 
algebraic, computation on the primitives. 

4.6 Handling boundary conditions 

MA treats the boundary conditions in a generic way. At 
first it identifies the type of boundary condition. In case of 
cauchy conditions the equation is discretized and simplified 
after substituting the grid parameters ; discretized equation for 
internal grid points is equated with the discretized boundary 
condition for eleminating the fictitous grid point and the final 
boundary equation is sent to the C module through a recursive 
data structure. In case of elliptic problenis special measures 


45 



bl-ioiild be t«t|.en for the corner nodes. Corvier nodes sharing 
neumaum conditions a^re equaited to e'J iininate the fictitous grid 

p 0 i n 't h t» 


4.7 Explanation 

One of the importamt features of expert systems is 
tho capability of answering the user query. It has both 
psychological avnd rnateria'l importance; it often give the user a 
feeling of assistance not as of being commanded by the systeni. 
MoreciVf-r system c'xplanationn are useful for a novice, who should 
be able to look through the system. 

The t,omp]et(? systeru is classified into several levels. A track of 
all the lull's fired at every level is maintained. When the user 
type why, system gets the recent rule and level from the database 
and answers it ’intelligently'. Although the technique for such 
quest :i ore aiiswe r ing is not very sophisticated, for our purpose it 
does the job effectively, as the user session is comparatively 
siirtpler here. 
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Chaftei' fjvt' Concluiiiorri 

5.1 Discussion 

As a v<.*ry limited number of numerical ei-'.i;>oi't syoteriis arc 
availtUjjc for persona] computer usei s MA ran bo of much use in 
the scientific and engineering af.'pl i cat ions , MA can handle a wide 
range of FDf"o with different kinds of boundary conditions. 
Presently available state- of-ar t analysis softwares CANSYS, 
NASTPAN, IDF~A9 etc. 3 can iiandle only specific types of PDPs. MA 
tertainiy overcomes that weakness by iuidressing the problem in 
-I generii: w<iy. Tlte knowlwdge acquisition module of MA makes it 
dynamic to let a numerical analyst experiment with new schemes. 
The unique feature of handling directly the continuous physics 
problem is handy and approaches towards a natural way of 
thinking. I'he present prototype version of MA has several 
limitations too and those will be discussed in the following 
s e c 1 1 u 11 . 


5.2 Limitations and Scope for future work 

As mentiuneil earlier, the major limitation of MA is that it can 
only deal with rectangular domain. Provisions are kept for 
interfacing MA with a solid rnodula which will expand the problem 
domain. Some of the features of symbolic computation module are 
ad hoc? for example, linearizing of the nonlinear source term may 
soiuetime genercite even simplified expressions in a crude form. 
This is because it cannot put together complex similar terms. But 
this does not result errors in computation. So far the system is 
tested for a variety of problems. Graphics module needs more 
attention. At present MA is equiped with two dimensional 



gi'rtf..!.i. c‘-. . 31) graphicu ii, essential especially for plotting 

solutj.oii‘3 of elliptic equations. Error handling file is 
iru oiiipliJtf ; it needs a rigorous long system testing along with 
the user f (Muibacks . 

Conclusion 

A new problem solving environment, MA for solving PDEs is 
proposed and implemented. This work several unique features along 
witii sumo drawbacks, discussed earlier. Earlier works on PDE 
solver were limited eitlier by their too rigid handling of the 
problerii (procedural subroutines) or by their inefficient nuinber 
crunching power (systeiiis written completely in relational 
langut^ges ) . Although some expert systems for solving PDEs, like 
AGNEP), riDISOL, Elliptic Expert , PLOD/PC etc. are capable of 

tacf.ling a wide range of PDEs, but all of them generate 
pra< odural codes, winch need to carry a procedural compiler 
(either internally or externsilly) and thereby demand more disk 
and inemory space; moreover this is time consuming too- MA has 
overcome tfris difficulty. All oF the above mentioned systems 
addres'a tlie pcobleiit as a PDE not as a phenomena. This is really 
an unwantcui feature of an expert system, which is supposed to 
provide a higher level of user interface. MA certainly got rid of 
this limitation. Although the storage of continuous physics 
knowledge is somewhat ad hoc, this can be overcome by deriving 
governing PDEs from the fundamental equations using symbolic 
operators . 
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Fie. 6 feedback of MA -preprocessor 


S4 


Welcome to EXPERT SYSTEM MA 
This is an expert system for solving 
Partial Differential Equations 

All copyrights reserved to 
Hillol Kargupta 
IIT Kanpur, India 


Fig.7 - System plate no.1 
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Fig. 1l - System plate no. 4 
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Please, define the unknown ternn(s) 
found in the equation. Note that, 
unknown terms can either be a 
function of the dependent variable 
and independent variables or constants. 


diffusivity= 0.01 


|Esc - Mam mentil j? - Help | 


laf »Tmi 


|why - Ext.lairi 


Fig. 12 - System plate no. 5 
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EQUATION TO BE SOLVED : 0.01*d2u/dx2-du/dt=0 

THIS IS A LINEAR, PARABOLIC IN SPACE, FIRST ORDER IN TIME, 
AND SECOND ORDER IN SPATIAL COORDINATES 

We need TWO boundary conditions along the X-axis and 
ONE initial condition to solve the problem. 

I will be asking you about the boundary & initial conditions 
and the grid parameters, soon. 

Press any key to proceed 


E sc - Main menu | |? - HelfT] - QuiFI jwhy - Explain j 


Fig. 13 - System plate no. G 
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Enter X direction width =1 

Enter total number of X grid points, Nx =11 

Enter time step, Dt =0.001 

Enter total number of time steps, Nt =30 

Enter initial distribution, Uo =100 

First boundary condition =1000 

Second boundary condition =100 


irs~l^in mg77i^ |? - Herp [ 


tamirrti 


11 Why - Explain j 


Fig. 14 - System plate no. 7 
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Fig. 15 - System plate no. S 
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6.1 Getting Started s 

Tills chapter describes the basic operation of the expert system 
MA, including how to use the menu systein and control the 
env i r onment . 

Tite README! file on the distribution disk gives a complete list of 
tl'ie files supplied on the distribution disks. MA comes ready 
install and run on an IBM PC or fully compatible computer. 

F’ollow one of these two procedures to load MA into your system : 

# If you are running MA from floppy disks s 

1. Turn on and boot up your computer. 

2. Take your system start-up disk out of the dist 
drive. 

3. At the A!\>prompt, type ma and press Enter. 

# If you are running MA from your hard disk s 

1. Turn on your computer. 

2. Log onto the hard disk and change to the directory 
where MA resides. 

3. Type ma and press Enter. 

You should now see the logon message shown in Fig. 7. Press any 
character key; MA’s inain window and menu will appear, as shown in 
Fig. 8. 

6.2 A Quick Guide to the Menus and Hot Keys 

In this section an introduction to MA's menus and some special 
shortcut keys (known in the vernacular as hot keys) is given. 
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Tho Main Menu : 


MA ' y fitain menu showy you the commands and pulldown menus 
available. You can select an item from tlie main menu in one of 
the Following two ways s 

1. When main menu is active, by pressing the associated 
hiyi'il i.gh ted capital letter (S for Solve, 0 for Opt,ions, 
etc . ) 

2. When main menu is active, by first moving the highlight 
bar with the arrow keys, and then pressing Enter 


The Pull-down Menus! 

Four of the items on the mam menu invoke pull-down menus when 
you press them? these pull-down menu items are Solve, Utilities, 
Database and Options. Here’s a brief description of the main menu 
items : 

a) Solve 


Using this item you can invoke the problem solving 
riiodule, which initiates a course oF dialog with the user 
and finally solves the problem. Two options are of-Pered, 
natiiely ! 

Feed equation - You can directly feed the PDE to 
be solved using this facility. 

State problem - The physical problem can bpstated , 
which in turn fetches the PDE and relevant informations 
from the knowledge base. 



b) 


Utilities 


The items on this menu allow you to handle •Tiles, 
manipulate di recto ry, invoke DOS, call graphics modiile. 
Ihe listing o-T ti'ie items are s 

1. Copy file One can copy a file to a destination file by 

using this option. 

?. Print file This can be invoked to get the printout of a 

file directly fiom tho environment. 

3. OS shell One can go to tlie operatiivg system directly 

from the system? type exit to come back 
to the systeni. 

4. Graph Plot One tan directly call the graphics module 

for plotting graphs, by stating tlie 
graphics filename. 

5. Sketch Pad At present no interface is provided 

with MA. This is kept for future 

modifications, which will enable the 
user to define the problem domain 

graphi cal ly . 

c) Database 

Database option handles the manipulation of the 
knowledgebase. Four choices are available: 

1. Display database file - All the updatable knowledge 
is stored in tlie file mabase.dba. This option allows the 


7 1 



to J. not' at tho without allowing oriy 

change. Press Esc key for conring out from display inode. 

2. Update database ~ System t nowl edgebass can bo 
uijdcited using this option. At first system will ask you 
about tlie physical name of the problem, eg. heat 
conduction, vibration etc. and then the governing PDE. 
Var iable properties can be introduced. Use default 
denotations for variables, otherwise define new 
variables using menu Options/Variable denotations. 
Remember once you update the database all the new 
informations are kept in RAM only, not in the file. So 
you must invoke Save database to file menu if you want 
to retain it. 

3. Edit database file - Any bug in the database can be 
eleminated using this option. This option invoke an 
editor loaded with the database file "mabase . dba" . 

4. Save database in file ~ This option saves the 
current database in the file "mbase . dba" . It always 
creates a back-up, named mabase. bak. 

d) Options 

This menu provides several subsequent pull-down menus, 
from which you can change the default description of the 
dependent and indepentent variables, computational 
scheme and discretization scheme. 



■1 . Variable denotation 


Dif-feri^nt Uemutat 3 oni^ for differonl 


type of problenis makes the system handy 
and PleKible. Default denotation for the 
dependent variable is u, K-a^is x, and 
y-axls y. 

2. Discretization scheme iiser can specify the discretization 

scheme beforehand. This can also be dune 
duriny the dialogue session with the 
systein. At present Discretization scheme 
offers a choice among the Fully implicit 
scheme, Crank-Ni colson scheme and 
explicit scheme. 

3. Computational scheme Similarly the computational scheme 

can also be beforehand. At present 
Succesive Over Relaxation (SOR), Gauss- 
Siedel, Tri-diagonal algorithm, and 
Gauss-Jordan methods sire available for 
solving the set of algebraic equations. 


e) Quit 

This item should be chosen to come out of 

the system. 


Hot Keys and commands 

esc go to main menu 

? invoke on line help file 

why query about the current process 

quit 
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6.3 Solving a Partial Differential Equation s 


Solve ofitioTi provides two choices fur solving the 
pi'oblofii. The State problem option invokes the 
preprocessor directly and displays a pop-~up menu, 
containing all the physical problems, defined in t!'ie 
database, and then asts the user about the unknowns 
present in the equation, if any. Fig. 9 shows the 
procedural flow of MA. 

If the Feed equation option of the main menu is hit, systeiit 
will want the equation to be solved, by prompting "Enter 
the equation to be solved", as shown in Fig 10; after 
piarsing the equation the list of unidentified terms are 
sent back to the user for proper defination. Once the 
equation is completely defined the processed equation is 
sent back to the user for approval along with several 
useful informations about it. This is shown in Fig. 6. 
If the discretization scheme is not chosen beforehand a 
menu offering several discretization schemes appears in 
the screen. After this the dialog session for obtaining 
the grid parameters and boundary conditions comes into 
action. After the boundary conditions are identified and 
processed, discretized coefficients are sent back to the 
user. Then system suggests the best computational 
scheme; but the user can deny' it and choose the desired 
one from the menu, as shown in Fig 6. At any stage 
user may query the system by typing why. Now the control 
goes to the procedural number crunching modules. Output 
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re<iu3t is stored in both an output file and a graphics 
data -File, Results are also displayed on the screen at 
the sarne time. User can stop the number crunching 
process by pressing Ctrl— break key. As soon as the 
computation work is over message window prompt will let 
it know. A pop-up menu af>pears in the mid-screen which 
offers several actions- Options are equation dependent. 
Followings are the different options available for 
different type of equations 5 

a. Elliptic equation 

1) U (dependent variable) vs. X (jt-axis variable) 
If this option is chosen the y-values for which a U/X 
graph will be plotted, are further askied. 

2) U (dependent variable) vs. Y (y-axis variable) 
If this option is chosen the 5{“values for which a U/Y 
grapl ‘1 will be plotted, are further asked. 

3) Dad- to main menu 

b. F’arabolic equation 

1) Sitnple plot - This option offers an usual grapdi 
plotting with the spatial domain along x axis and the 
dej>etident vaiiblo along y-axis. 

2) Back to Ilia in iiienu 
b. Hyperbolic equation 

1) Simple plot - litis option offers an usual graph 
plotting with the spatial domain along x axis and the 
dependent varible along y--axis. 


2) Back to main menu 



6.5 Mathemat i cal functions 


Several mathemat i cal functions are avail lable in 
the library. A listing of these functions are s 


s i n ( X > 

sin function 

c.05(X) 

coo function 

tan(X) 

tan function 

cose c ( X ) 

cosec function 

se c ( X ) 

sec function 

cot<X) 

cot function 

ln(X) 

natural logarithm of X(base e) 

logCX) 

base 10 logaritlim of X 

e X p ( X ) 

e raised to the value of X 

abs(X) 

absolute value of X 

X-'Y 

X to the power Y 


6.6 Errors 

A nuiiiber of error^may occur at run time. All the 

errors aiiong with tlieir reasons are stared in MA.ERR. In case of 
any error system sents a message to the message window. If the 

user commits a^ny mistake in typing either tl'ie aquation or the 

paramc'tors , system points out tlie mistake immediately and lets the 
user correct it. In case of inemory overflows, which may occur due 
to ei'.tensive grid r esol ut i on , ve ry large number of iterations etc., 
system sends a message in the message window and resets- the 

system. 





Appendix B ■“ Case study 


Expert system MA has a robust FDM solver, which can handle a 
wide ravige of partial differential equations. Our system has 
been tested for a wide range of PDFs and so far at has proved 
its computational ability. In this section we shall discuss a 
sample problem, to illustrate the complete process of solving a 
PDF using MA. Although a simple problem is taken here, MA can 
handle complex situations too. 

Consider the following heat transfer probloni s 

a) The wall of a furnace is 1ft thict and is made of brick, which 

has a thermal diffusivity of 0.01 ft'^/h. The temperature of the 
wall 35 100 °F when the furnace is off. When the furnace is 

fired, the temperature on the inside face of tl'ie wall reaches 
lOOO^r quite rapidly. The temperature of tl'ie outside face of the 
wall IS maintained at 100*^F by natural convection. Determine the 
evolut'ion of temperature profiles within the brick wall. 

b) Insulatioji is placed on the outside surface of tl'ie wall. 

Assume this is perfect insulation and show the evolution of the 
temperature profiles within the wall when the furnace is fired to 
lOOO^F. 

Both the above problems can be solved by MA 
with ease. F'ig.11 shows tl'ie problem type menu, whicl'i comes as 
soon as one hits the State problem option. Heat conduction (ID) 
can be chosen, which results in the preprocessing session, as 
shown in Fig. 12, Tht' System then asks for the value of diffusivity. 
Once all the terms of the PDE, governing the problem are defined 
a feedback message is sent to the user. This is shown in Fig. 13. 
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Now the I’lrocessor module is invoked, which asks the user about 
the yrid parameters and boundary conditions ai shown in Fig.’!-*!. 
All the values are entered as stated in tiie -first probleni? 
finally the computation module is invoked, whicl'i saves the output 
to the user specified file, displays on the screen and also keeps 
in a database cliain for graphics. Graphics output of the first 
problem is shown in Fig. 16. 

Fig. 15, Fig. 17 show the changed boundary condition and the 
CO r res ponding solution in graphics, respectively. 



